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SCIENTIFIC  ACCOMPLISHMENTS 


The  goal  of  this  program  has  been  to  develop  theoretical  understanding  of  processes  of 
decoherence  and  relaxation,  as  well  as  of  qubit  interactions,  in  order  to  evaluate  quantum 
computing  designs  in  solid  state,  with  emphasis  on  spin-based  schemes  in  low- 
temperature,  high  magnetic  field,  quantum  Hall  effect  regime  semiconductor 
heterostructures.  A  broader  aim  has  been  to  contribute  to  the  understanding  of  spin 
control  and  measurement  in  semiconductor  and  other  materials,  to  support  experimental 
efforts  in  quantum  information  processing  and  spintronic  device  development. 

For  a  quantum-computing  design  based  on  quantum-Hall  semiconductor  heterostructures, 
we  have  accomplished,  in  several  publications,  calculation  of  all  the  relevant  control  and 
relaxation/decoherence  time  scales.  This  has  allowed  us  to  combine  elements  of  the  1998 
quantum  computing  designs  by  Privman,  Vagner  and  Kventsel  [1],  and  by  Kane  [2],  with 
the  new  idea  [3]  of  nuclear-spin  qubit  interactions  mediated  indirectly  via  the  bound  outer 
electrons  of  impurity  atoms  whose  nuclear  spins  1/2  are  the  qubits.  These  electrons,  in 
turn,  interact  via  the  two-dimensional  electron  gas  in  the  quantum  Hall  effect  regime.  The 
resulting  quantum  computing  scheme  [3-5]  retains  all  the  gate-control  and  measurement 
aspects  of  the  proposal  by  Kane,  but  allows  qubit  spacing  at  distances  of  order  100  nm, 
attainable  with  the  present-day  semiconductor-heterostructure  device  technologies.  Two 
articles  are  attached  detailing  these  results:  Appendix  A  is  an  overview-style  paper  [5], 
whereas  Appendix  B  is  a  more  technical  work  [3]. 

In  order  to  generally  gain  understanding  of  decoherence  in  materials  considered  for 
quantum  computing  realizations,  we  [6],  as  well  as  other  groups  [7,8],  have  studied  the 
low-temperature  dynamics  of  a  shallow  donor  P  impurity  or  quantum-dot  bound- 
electron  spins  in  silicon.  Specifically,  we  considered  its  interactions  with  the  bath  of 
nuclear  spins  of  the  29Si  isotope.  For  small  applied  magnetic  fields,  the  electron  spin 
relaxation  is  controlled  by  the  steady  state  distribution  of  the  nuclear  spins.  We  calculated 
the  relaxation  times  7]  and  T2  as  functions  of  the  external  magnetic  field,  and  concluded 
that  nuclear  spins  play  an  important  role  in  the  donor  electron  spin  decoherence  in  Si:P  at 
low  magnetic  fields. 

For  general  evaluation  of  relaxation  processes  at  short  times,  of  relevance  for  gate  control 
of  qubits  and  qubit-qubit  interactions,  we  have  initiated  development  of  a  new 
approximation  scheme  [9,10]  for  evaluation  of  decoherence.  At  low  temperatures,  the 
approximation  is  argued  to  apply  at  intennediate  times  as  well,  up  to  the  thermal  time 
scale  fil  kT .  It  then  provides  an  approach  complementary  to  Markovian-type 
approximations,  and  is  appropriate  for  evaluation  of  deviations  from  pure  states  in 
quantum  computing  models.  Our  work  to  apply  this  method  to  qubits  in  quantum¬ 
computing  architectures,  is  ongoing,  continuing  under  the  new  ARO  grant.  In  various 
collaborations,  we  have  also  contributed  to  studies  of  spin  transport,  with  applications  in 
spintronics  and  quantum  measurement:  the  full  list  of  publications  is  provided  in  a  later 
section.  Appendix  C  is  a  manuscript  [10]  detailing  our  new  approach  to  evaluation  of 
relaxation  and  decoherence  at  short  times. 
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Abstract 

We  survey  recent  work  on  designing  and  evaluating  quantum  computing  implementations  based  on  nuclear  or  bound-electron 
spins  in  semiconductor  heterostructures  at  low  temperatures  and  in  high  magnetic  fields.  General  overview  is  followed  by  a 
summary  of  results  of  our  theoretical  calculations  of  decoherence  time  scales  and  spin-spin  interactions.  The  latter  were  carried 
out  for  systems  for  which  the  two-dimensional  electron  gas  provides  the  dominant  carrier  for  spin  dynamics  via  exchange  of 
spin-excitons  in  the  integer  quantum  Hall  regime.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 

PACS:  73.20.Dx;  71.70.Ej;  03.67.Lx;  76.60.-k 
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1.  Introduction 

The  field  of  quantum  computing  has  seen  explo¬ 
sive  growth  of  experimental  and  theoretical  interest. 
The  promise  of  quantum  computing  [1-5]  has  been  in 
exponential  speedup  of  certain  calculations  via  quan¬ 
tum  parallelism.  In  Fig.  1,  the  top  flow  chart  shows 
the  “classical”  computation  which  starts  from  binary 
input  states  and  results  in  binary  output  states.  The  ac¬ 
tual  dynamics  is  not  really  that  of  Newtonian  classi¬ 
cal  mechanics.  Rather  the  computation  involves  many- 
body  irreversible  “gate”  device  components,  made  of 
semiconductor  materials  in  modern  computers,  which 
evolve  irreversibly,  “thermodynamically”  according  to 
the  laws  of  statistical  mechanics.  As  the  size  of  the 
modern  computer  components  approaches  atomic,  the 
many-body  quantum  behavior  will  have  to  be  ac¬ 
counted  for  in  any  case  [6]. 
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The  idea  of  quantum  computing,  however,  is  not 
just  to  account  for,  but  to  actually  utilize  the  quantum- 
mechanical  dynamical  behavior.  This  is  not  an  easy 
task.  Quantum  mechanics  allows  for  parallelism  in 
evolution:  one  can  “process”  a  linear  superposition 
of  several  input  states  at  once,  as  illustrated  in  the 
lower  flow  chart  in  Fig.  1.  The  price  paid  is  that 
coherent  processing  of  information,  according  to  the 
law  of  quantum  mechanics,  must  be  accomplished 
in  systems  much  larger  than  atomic-size  (or  more 
importantly,  with  many  degrees  of  freedom).  There 
are  numerous  conceptual  and  experimental  obstacles 
to  accomplishing  this  task,  that  have  generated  a  lot 
of  interest,  excitement,  and  new  results  in  computer 
science,  physics,  and  engineering. 

The  functioning  of  a  quantum  computer  involves 
initialization  of  the  input  state,  then  the  actual  dy¬ 
namical  evolution  corresponding  to  computation,  and 
finally  reading  off  the  result.  Various  specific  re¬ 
quirements  for  implementation  have  been  identified 
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Fig.  1 .  Comparison  of  the  classical  and  quantum  approaches  to  com¬ 
puting.  The  upper  flow  chart  schematically  represents  implemen¬ 
tation  of  a  traditional  irreversible  “classical”  computation  process, 
where  transformation  of  the  input  set  of  bits  into  the  result  is  ac¬ 
complished  by  a  succession  of  irreversible  gates.  Owing  to  their  irre¬ 
versibility,  the  gates  can  be  connected  in  space  rather  than  switched 
on  and  off  at  different  times.  The  lower  flow  chart  shows  quantum 
processing  of  information,  where  the  input  and  the  final  result  are 
both  in  superposition  states,  yielding  quantum  parallelism.  The  dy¬ 
namics  is  reversible:  there  is  a  one-to-one  correspondence  between 
the  initial  and  final  states.  Therefore,  number  of  the  input  and  out¬ 
put  quantum  bit  (qubits)  is  the  same  even  though  some  of  the  output 
qubits  (set  in  a  smaller  font)  might  not  be  used  in  the  final  extrac¬ 
tion  of  the  classical  result  by  measurement.  The  quantum  gates  are 
applied  in  succession  by  being  switched  on  and  of  at  different  times 
during  the  computation.  The  question  mark  signifies  the  difficulty  of 
finding  quantum  algorithms  that  retain  the  power  of  quantum  paral¬ 
lelism  after  measurement  needed  to  read  off  the  final  result  as  clas¬ 
sical  information. 


[2-5];  here  we  provide  only  a  limited  introductory 
overview. 

Let  us  begin  by  considering  the  reading  off  of  the 
final  result.  The  reason  for  the  question  mark  in  the 
lower  chart  in  Fig.  1  is  that  quantum  measurement  of 
the  final  superposition  state  can  erase  the  gain  of  the 
parallel  dynamics,  by  collapsing  the  wave  function. 
Therefore,  a  key  issue  in  quantum  computing  has 
been  to  find  those  algorithms  for  which  the  readout 
of  the  final  state,  by  way  of  projecting  out  a  certain 
average  property,  still  retains  the  power  of  the  quantum 
parallelism.  To  date,  only  few  such  examples  are 
known  [1, 3,4,7],  the  most  celebrated  being  the  Shor 
algorithm  [1]  for  factoring  of  integers,  the  invention  of 
which  boosted  quantum  computing  from  an  obscure 
theoretical  field  to  a  mainstream  research  topic. 

The  preparation  of  the  initial  state  does  not  seem 
to  present  a  problem  for  most  quantum  computing 
realizations  [2-5],  except  perhaps  the  ensemble  liq¬ 


uid  state  NMR  approach  [8,9]  which  relies  on  the 
initial  thermal  distribution  to  produce  deviation  of 
the  density  matrix  from  the  equal-probability  mixture 
state.  In  most  other  approaches,  the  initial  state  can 
be  produced  by  first  fully  polarizing  the  quantum  bits 
(qubits),  i.e.  putting  them  in  one  of  the  two  quantum 
levels.  Note  that  we  consider  two-state  qubits  here,  re¬ 
alized,  for  instance,  by  spins  1/2  of  nuclei  or  gate- 
or  impurity-bound  electrons,  in  applied  magnetic  field. 
The  fully  polarized  state  is  then  subject  to  gate  opera¬ 
tions  to  form  the  desired  input  state.  Part  of  a  quantum¬ 
computing  algorithm  should  be  the  prescription  on 
how  to  choose  the  initial  state  to  represent  the  clas¬ 
sical  information  of  the  input,  like  the  input  integer  in 
the  factoring.  In  most  cases,  this  prescription  is  easily 
accomplished  by  single-qubit  and  two-qubit  gates. 

The  actual  dynamical  evolution  (the  process  of 
computation)  in  quantum  computing  is  fully  reversible 
and  nondissipative,  unlike  classical  computing.  Much 
progress  has  been  made  in  resolving  both  the  con¬ 
ceptual  and  computer-engineering  “design”  issues  for 
quantum  computation.  Specifically,  the  computation 
can  be  carried  out  [2-5,10-13]  by  a  universal  set  of 
gates;  single-qubit  rotations  and  nearly  any  two-qubit 
gate.  The  gates  are  not  connected  in  space  like  in  clas¬ 
sical  computers  but  are  activated  in  succession  in  time, 
to  control  single-spin  dynamics  and  also  switch  on  and 
off  two-spin  interactions  (we  use  “spin”  and  “qubit” 
interchangeably). 

Many  interesting  matters  have  been  resolved,  which 
are  not  reviewed  here.  These  include  the  understand¬ 
ing  of  how  the  finiteness  of  the  state  space  (i.e.  two 
states  for  spin  1  /2)  replaces  the  “classical”  digitaliza¬ 
tion  in  quantum  computing.  Also,  the  “classical”  copy¬ 
ing  (fan-out)  function  is  not  possible  in  quantum  me¬ 
chanics.  It  is  replaced  by  entanglement  with  ancillary 
qubits  to  accomplish  redundancy  needed  for  error  cor¬ 
rection  [14-20],  Sources  of  errors  due  to  interactions 
with  environment  in  quantum  mechanics  involve  not 
only  the  usual  relaxation  (thermalization)  but  also  loss 
of  coherence  [21-28].  This  quantum  decoherence  (de¬ 
phasing)  can  be  faster  than  relaxation  because  it  does 
not  require  energy  exchange. 

A  conceptually  important  issue  has  been  the  scala¬ 
bility  of  quantum  computing:  can  one  process  macro¬ 
scopic  ally  large  amounts  of  information  by  utilizing 
quantum  error  correction  based  on  redundancy  via  en¬ 
tanglement  with  ancillary  qubits?  The  affirmative  an- 
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swer  to  this  question  has  been  one  of  the  triumphs  of 
the  theory  [14-20],  It  provided  a  new  paradigm  for 
emergence  of  controlled/organized  macroscopic  be¬ 
havior  from  microscopic  dynamics,  on  par  with  the 
conceptual  possibility  of  living  organisms,  which  we 
observe  but  cannot  yet  “manufacture”,  and  million- 
gate  classical  computers  which  are  man-made. 

With  all  these  theoretical  advances  at  hand,  the  next 
step  is  to  ask  whether  a  man-made  quantum  computer 
can  be  realized?  There  have  been  several  experimental 
directions  of  exploration,  most  presently  are  still  at  the 
level  of  one  or  two  qubits,  or,  for  ensemble  liquid- 
state  NMR,  which  emulates  quantum  dynamics  by 
evolution  of  the  density  matrix  of  a  large  collection 
of  molecules,  5-7  qubits. 

In  this  introductory  survey,  we  summarize  results 
of  our  work  on  two-spin  interactions  and  spin  decoher¬ 
ence  in  semiconductor  heterostructures.  In  Section  2, 
we  consider  the  spin-based  quantum  computing  pro¬ 
posals  in  such  systems.  Time  scales  of  relaxation  and 
decoherence  are  addressed  in  Section  3.  Finally,  Sec¬ 
tion  4  reports  results  for  models  with  nuclear  spins  as 
qubits. 

2.  Spin-based  quantum  computing  in 
semiconductor  heterostructures 

The  general  layout  of  a  solid-state  quantum  com¬ 
puter  is  shown  in  Fig.  2.  Qubits  are  positioned  with 
precision  of  few  nanometers  in  a  heterostructure.  One 
must  propose  how  to  effect  and  control  single-qubit 
interactions,  two-qubit  interactions,  and  explore  how 
the  controlled  dynamics  owing  to  these  interactions 
compares  to  decoherence  and  relaxation.  The  proposal 
must  include  ideas  for  implementation  of  initializa¬ 
tion,  readout,  and  gate  functions. 

The  first  proposal  including  all  these  components 
was  for  qubits  realized  in  an  array  of  quantum  dots 
[29]  coupled  by  electron  tunneling.  The  first  spin- 
based  proposal  [30]  utilized  nuclear  spins  coupled  by 
the  two-dimensional  electron  gas,  the  latter  in  the  dis¬ 
sipationless  integer  quantum  Hall  state  [31]  that  re¬ 
quires  low  temperatures  and  high  magnetic  fields.  An 
important  advancement  was  the  work  of  Kane  [32] 
where  gate  control  of  nuclear-spins  of  donor  impuri¬ 
ties,  separated  less  than  10  nm  and  coupled  via  the 
outer  impurity  electrons  which  are  bound  at  low  tem- 
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Fig.  2.  Schematic  illustration  of  a  semiconductor  heterostructure 
quantum  information  processor.  The  qubits,  represented  by  the 
arrows  overlaying  heavy  dots,  are  spins  1  /2  of  nuclei  or  localized 
electrons.  Individual  control  of  the  temporal  evolution  of  the  spins 
can  be  achieved  with  the  use  of  external  electromagnetic  radiation, 
i.e.  NMR  or  ESR  pulses.  The  spins  are  also  coupled  with  each 
other  via  interaction  mediated  by  the  two-dimensional  electron 
gas  in  the  heterostructure,  or  by  other  means.  The  external  and 
internal  interactions  can  be  controlled  by  gates  formed  on  top  of 
the  heterostructure.  The  external  environment,  that  includes  crystal 
lattice,  electron  gas,  defects,  impurity  potentials,  causes  relaxation 
and  decoherence  of  the  qubits. 

peratures,  was  proposed.  Most  of  these  ideas  also 
apply  to  electron-spin  qubits,  bound  at  impurities, 
in  quantum  dots,  or  directly  by  gates.  Several  elab¬ 
orate  solid-state  heterostructure  quantum  computing 
schemes  have  been  proposed  in  the  literature  recently 
[28,33-41].  There  are  also  other  promising  proposals 
involving  surface  geometries:  superconducting  elec¬ 
tronics  [42 — 46]  and  electrons  on  the  surface  of  liquid 
helium  [47], 

There  have  been  several  planned  and  ongoing 
experimental  efforts  [32-36,43^-5,48-54]  ultimately 
aimed  at  solid-state  quantum  computing  and  other 
quantum  information  processing  realizations.  The  fi¬ 
nal  geometry  is  expected  to  be  most  sensitive  to  the 
implementation  of  readout,  because  it  involves  quan¬ 
tum  measurement,  i.e.  supposedly  interaction  with 
or  transfer  of  information  to  a  macroscopic  device. 
Therefore,  much  of  the  experimental  effort  presently 
has  been  focused  on  single-qubit  (single-spin)  mea¬ 
surement  approaches. 

The  theoretical  efforts  can  be  divided  into  two  ma¬ 
jors  tasks.  The  process  of  single-spin  measurement 
must  be  understood  for  the  readout  stage  of  quantum 
computing.  Several  conceptual  and  calculational  ad¬ 
vancements  have  been  made  in  understanding  quan¬ 
tum  measurement  [26,32-36,46,50,55,56]  as  it  applies 
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to  atomic-size  qubit  systems  interacting  with  environ¬ 
ment  and  typically  “measured”  directly  by  the  effect 
of  the  spin-qubit  state  on  transport,  or  first  transferring 
the  spin  state  to  a  charge  state  that  is  easier  to  measure, 
e.g.,  in  single-electron  transistors  and  similar  devices. 

In  this  survey,  we  outline  results  of  the  second 
evaluation  task:  that  of  understanding  the  processes 
and  times  scales  involved  in  the  dynamics  of  the  actual 
computation.  As  summarized  in  Fig.  3,  this  main 
stage  of  the  quantum  computation  process  involves 
control  of  spins  and  their  interactions.  It  also  involves 
processes  that  we  do  not  control  and  are  trying  to 
minimize:  relaxation  and  decoherence. 

Control  of  individual  qubits  is  usually  accom¬ 
plished  externally.  For  nuclear  spins,  NMR  radio¬ 
frequency  radiation  can  be  used,  see  Fig.  2.  For  elec¬ 
tron  spins,  the  ESR  microwave  frequencies  are  suit¬ 
able.  Such  radiation  cannot  be  focused  on  the  scale 
of  10-100  nm.  Instead,  selectivity  must  be  accom¬ 
plished  by  independent  means.  Several  proposals  ex¬ 
ist,  the  most  promising  being  control  by  gates.  The  ap¬ 
plied  gate  voltage  modifies  the  electronic  wave  func¬ 
tion  changing  interactions  and  therefore  resonant  fre¬ 
quencies.  We  will  denote  the  time  scale  of  the  external 
single-qubit  control  by  Text.  This  can  be  the  Rabi  time 
of  a  spin  flip. 

The  qubit-qubit  interactions  are  typically  assumed 
to  be  mediated  by  electrons  that  “visit”  both  qubit 
environments.  For  instance,  in  liquid-state  ensemble 
NMR  [8,9]  with  complex  molecules,  or  in  the  original 
model  [32]  of  phosphorous  impurity  donors  in  silicon, 
the  wave  functions  of  the  valence,  outer  electrons  of 
nearby  qubits  overlap.  Specifically,  in  the  P  donor 
case,  the  single  outer  electron  of  the  donor  atom 
remains  bound  at  low  temperatures  but  has  orbital 
radius  of  order  2  nm  owing  to  the  large  dielectric 
constant  of  the  silicon  host.  Therefore,  it  is  hoped 
that  these  electrons,  in  nearby  donors  positioned  as  in 
Fig.  2,  will  mediate  nuclear-spin  qubit  interactions. 

Our  approach  [27,28]  allows  for  larger  qubit  sep¬ 
aration,  up  to  order  100  nm,  by  relying  on  the  two- 
dimensional  electron  gas  in  the  heterostructure  to  me¬ 
diate  qubit-qubit  interactions.  This  two-dimensional 
electron  gas  is  usually  obtained  by  spontaneous  or 
gate-induced  transfer  of  electrons  from  impurities 
to  the  two-dimensional  interface  layer  in  which  the 
qubits  are  positioned.  The  source  impurities  are  lo¬ 
cated  at  some  separation  from  this  layer  or  in  the 


bulk.  The  two-dimensional  electron  gas  can  be  made 
nondissipative  in  certain  ranges  of  large  applied  mag¬ 
netic  fields  at  low  temperatures,  when  these  conduc¬ 
tion  electrons  in  the  layer  are  in  the  integer  quantum 
Hall  effect  state.  Owing  to  this  property  and  also  larger 
qubit  separation  allowed,  we  consider  this  the  most 
promising  approach  and  focus  our  present  review  on 
such  systems. 

The  time  scale  of  the  qubit-qubit  interactions  will 
be  denoted  by  7,m.  This  is  the  time  it  takes  to  accom¬ 
plish  a  two-qubit  quantum  gate,  such  as  CNOT  [2- 
5,57].  Typically  for  semiconductor  quantum  comput¬ 
ing  proposals,  7]nt  <  7’cxt,  and  in  fact  the  case  with 
Tint  Text  has  some  advantages  because  one  can  use 
several  fast  single-spin  flips  to  effectively  switch  in¬ 
teractions  of  some  qubits  off  over  the  gate  cycle.  An¬ 
other  approach  to  controlling  (on/off)  of  the  two-qubit 
interactions  is  by  gates,  see  Fig.  2,  which  affect  the 
two-dimensional  electron  gas  and  the  localized  elec¬ 
tron  wavefunctions. 

However,  the  same  conduction  electrons  that  pro¬ 
vide  the  qubit-qubit  interactions,  also  expose  the 
qubits  to  the  environment,  causing  relaxation  and  de¬ 
coherence.  Other  interactions  will  also  be  present,  that 
play  no  role  in  the  useful  quantum-computing  dynam¬ 
ics  but  contribute  to  these  undesirable  processes.  Re¬ 
laxation  and  decoherence,  and  their  associated  time 
scales,  are  addressed  in  the  next  section. 

3.  Time  scales  of  relaxation  and  decoherence 

The  processes  of  relaxation  and  decoherence  con¬ 
sidered  here  [21-28]  are  associated  with  the  dynamics 
of  a  small,  few-qubit  quantum  system  as  it  interacts 
with  the  environment.  Ultimately,  for  a  large,  multi¬ 
qubit  system,  many-body  quantum  chaos-like  behav¬ 
ior  must  also  be  accounted  for,  and  some  advances  in 
model  system  studies  have  been  reported  recently  [5, 
58].  Our  discussion  here  will  be  for  the  few-qubit  case 
mostly  because  it  allows  more  system-specific  investi¬ 
gations  for  actual  quantum-computing  proposals. 

Dynamical  processes  that  are  unwanted  in  quantum 
computing,  because  they  result  from  the  environmen¬ 
tal  influences  rather  than  from  the  controlled  radiation 
pulses  and  gate  potentials,  can  proceed  on  various  time 
scales.  In  fact,  it  is  not  guaranteed  that  processes  of 
various  types,  relaxation/thermalization  vs.  decoher- 
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ence/dephasing,  can  even  be  unambiguously  distinctly 
identified. 

At  low  temperatures,  it  is  generally  hoped  that  ther- 
malization,  which  requires  transfer  of  energy,  slows 
down.  If  the  fastest  such  processes  proceed  on  time 
scales  of  order  T\ ,  then  this  time  increases  at  low  tem¬ 
peratures  because  there  are  less  excitations  (phonons, 
electron  gas  modes,  etc.)  to  couple  the  small  quantum 
system  to  the  rest  of  the  solid-state  host  material. 

On  the  other  hand,  processes  that  do  not  require 
flow  of  energy  to  or  from  the  environment,  can  still 
effect  the  phase  of  the  quantum-superposition  ampli¬ 
tudes  and  cause  decoherence.  These  processes  can 
thus  proceed  faster,  on  the  time  scale  TV  While  these 
comments  seem  to  suggest  that  Ti  ^  7),  there  is  no 
obvious  reason  to  have  generally  73  7/  at  low  tem¬ 

peratures. 

However,  if  the  spectrum  of  the  dominant  excita¬ 
tions  mediating  the  qubit  coupling  (both  to  each  other 
and  to  the  host  material)  has  a  gap,  then  we  expect  that 
all  the  relaxation  and  decoherence  processes  will  be 
suppressed.  Furthermore,  the  suppression  of  the  relax¬ 
ation  will  be  exponential,  with  the  Boltzmann  factor 
for  that  energy  gap.  Then,  Ti_  <5IC  T\  will  be  satisfied 
but  also,  more  importantly,  the  actual  values  of  both 
time  scales  will  be  inordinately  large.  This  was  found, 
theoretically  and  experimentally,  to  be  the  case  for 
the  integer-quantum-Hall-state  two-dimensional  elec¬ 
tron  gas  as  mediator  of  the  localized-spin  (nuclear, 
electronic)  coupling  in  semiconductor  heterostructures 
[27,28,59-63]. 

It  is  important  to  emphasize  that  relaxation  and  de¬ 
coherence  are  really  many-body  properties  of  the  sys¬ 
tem  plus  environment.  Entanglement  with  the  environ¬ 
ment  owing  to  the  unwanted  couplings  results  in  the 
small  quantum  system  having  no  pure  wavefunction 
even  if  initially  it  was  prepared  in  a  pure  state.  Instead, 
it  can  be  described  by  a  statistical  mixture  represented 
by  a  density  matrix,  once  the  environment  is  traced 
over. 

This  reduced  density  matrix  of  the  system  is  ex¬ 
pected  to  evolve  to  the  thermal  one  at  large  times.  The 
approach  to  the  thermal  density  matrix,  which  is  di¬ 
agonal  in  the  system-energy  basis,  defines  the  time 
scale  T\.  If  the  temperature  is  low  enough,  then  there 
is  the  expectation,  see  [25,26]  and  references  therein, 
that  for  some  intermediate  time  scales,  of  order  T2, 
the  density  matrix  becomes  nearly-diagonal  in  a  basis 
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Control  qubits: 
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T\ 

&  decoherence: 

Ti 

(<A) 
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Control  and 
evolve  in  time 

Text  >  Tj„t  «  7  i  .  7] 


Measure 

Fig.  3.  Evaluation  of  quantum  computing  models.  One  of  the 
criteria  for  feasibility  of  quantum  computing  in  a  given  physical 
system  is  the  possibility  of  initialization  of  the  qubits  in  the 
desired  superposition  state.  Another  important  design  consideration 
is  control  of  qubit  states  and  of  their  interactions.  In  order  to 
implement  quantum  computing  effectively,  the  time  scales  for 
realization  of  single  and  two-qubit  logic  gates,  Text  and  7int> 
respectively,  should  be  several  orders  of  magnitude  smaller  than 
the  time  scales  of  relaxation  and  decoherence,  7)  and  73 .  The 
relationships  between  these  time  scales  are  further  explained  in  the 
text.  Finally,  efficient  and  reliable  measurement  of  the  output  state 
of  the  qubits  is  required  for  reading  off  the  result  of  the  computation 
and  presently  represents  a  formidable  experimental  challenge. 

which  is  determined  not  by  the  systems  Hamiltonian 
(energy),  but  by  the  interaction  operator  with  the  en¬ 
vironment.  This  latter  process  corresponds  to  loss  of 
quantum  coherence. 

As  emphasized  in  Fig.  3,  evaluation  of  a  quantum¬ 
computing  proposal  requires,  among  other  things,  es¬ 
tablishing  the  relation  Text,  7int  <5jC  73 ,  T\ .  Owing  to 
calculational  difficulties,  the  single-qubit  times  T\  2 
will  usually  be  used,  though,  as  mentioned  earlier, 
some  study  of  the  multi-qubit  “quantum  chaos”  effects 
may  be  required.  For  spin-qubit  quantum  computing 
in  semiconductor  heterostructures,  the  relation  is  typ¬ 
ically  rext  <33  7]nt  S3  T2  <33  7),  so  the  issue  is  usually 
how  small  is  the  quality  ratio  Q  =  7’mt/ Tn. 

The  required  value  of  Q ,  needed  for  fault- tolerant 
quantum  error  correction,  depends  on  the  physical 
model  of  error  sources  and  can  be  as  small  as  Q  = 
10“6-10“4,  See  [15,18-20],  or  as  large  as  Q  =  1/2, 
see  [64].  For  the  systems  of  interest  to  us  here,  spin 
qubits  in  semiconductor  structures,  the  value  of  Q  = 
10-5  is  a  reasonable  working  estimate.  Thus,  we  seek 
systems/conditions  with  Tmt/T2  ^  10-5. 

4.  Results  for  nuclear-spin  qubits 

In  this  section  we  outline  results  for  models  of 
quantum  computing  with  nuclear  spins  as  qubits,  and 
with  coupling  mediated  by  the  two-dimensional  elec- 
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tron  gas  in  the  integer  quantum  Hall  effect  state  [27, 
28,30].  In  strong  magnetic  fields,  the  spatial  states  of 
the  electrons  confined  in  the  two-dimensional  layer  in 
which  the  qubits  are  placed,  see  Fig.  2,  are  quantized 
by  the  field  to  resemble  free-space  Landau  levels.  The 
lattice  potential  and  the  impurities  actually  cause  for¬ 
mation  of  narrow  bands  instead  of  the  sharp  levels, 
separated  by  localized  states.  As  a  result,  for  ranges  of 
magnetic  field,  the  localized  states  fill  up  while  the  ex¬ 
tended  states  resemble  completely  filled  integer  num¬ 
ber  of  Landau  levels.  These  states  are  further  Zeeman 
split  owing  to  the  electron  spin.  At  low  temperatures, 
one  can  find  field  values  such  that  only  one  Zeeman 
sublevel  is  completely  filled  in  the  ground  state. 

The  electronic  state  in  such  systems,  that  show 
the  quantum  Hall  effect  [31]  in  conductivity,  are 
highly  correlated  and  nondissipative.  If  nuclear  spins 
are  used  as  qubits,  i.e.  atoms  with  nuclear  spin  1/2 
are  sparsely  positioned  in  the  zero-nuclear  spin  host, 
such  as  the  zero-nuclear-spin  isotope  28  of  Si,  which 
constitutes  92%  of  natural  silicone,  then  their  zero- 
temperature  relaxation  will  be  significantly  slowed 
down:  experimentally,  T\  —  103  sec  [62]. 

Localized  spins,  both  nuclear  and  electronic,  inter¬ 
act  by  exchanges  of  spin  excitons — spin  waves  con¬ 
sisting  of  a  superposition  of  bound  electron-hole  pair 
states.  The  spectrum  of  these  excitations  [65,66],  ob¬ 
served  experimentally  in  [67],  has  a  gap  correspond¬ 
ing  to  the  Zeeman  splitting.  This  gap  is  the  cause  of 
slow  relaxation  and  decoherence.  The  exchange  of  vir¬ 
tual  spin  excitons  mediates  the  qubit-qubit  interaction 
and  also,  via  scattering  of  virtual  excitons  from  im¬ 
purity  potentials,  relaxation  and  decoherence  of  single 
qubits. 

The  original  proposal  to  use  nuclear  spin  qubits  di¬ 
rectly  coupled  by  the  two-dimensional  electron  gas 
[30],  required  positioning  the  qubits  at  distances  com¬ 
parable  to  several  magnetic  lengths.  The  latter  is  of  or¬ 
der  10  nm  for  magnetic  fields  of  several  Tesla.  The 
qubit-qubit  interaction  decays  exponentially  on  this 
length  scale.  Recently,  we  proposed  a  new  improved 
model  [28]  in  which  the  qubit  interactions  are  me¬ 
diated  via  coupling  of  the  two-dimensional  electron 
gas  to  the  outer  impurity  electrons.  This  applies  if  the 
atoms,  whose  nuclear  spins  are  the  qubits,  are  single¬ 
electron  donors  such  as  the  isotope  31  of  R  These 
phosphorous  impurities  were  originally  utilized  in  the 
model  of  Kane  [32]  where  they  must  be  actually  posi¬ 


tioned  at  separations  of  about  4  nm  for  the  wavefunc- 
tions  of  the  outer  electrons,  which  are  bound  at  low 
temperatures,  to  overlap  significantly. 

In  our  new  improved  model  [28],  with  nuclear 
spins  coupling  to  the  outer  bound  electrons  which,  in 
turn,  interact  via  the  two-dimensional  electron  gas,  the 
interaction  turned  out  to  be  of  a  much  longer  range  as 
compared  to  the  model  of  [32]:  the  qubit  separation 
can  be  of  order  100  nm.  Another  advantage  is  that 
gate  control  of  the  individual  qubits  and  of  qubit- 
qubit  interactions  is  possible.  We  have  carried  out 
extensive  perturbative  many-body  calculations  [27,28, 
30,68]  allowing  estimation  of  7]nt  and  7/  for  both 
the  original  quantum-computing  proposal  [30]  and  its 
improved  version  [28],  where  the  main  improvement 
is  in  the  possibility  of  the  gate  control  along  the 
lines  of  [32],  The  “clock  speed”  of  the  improved 
model  is  also  faster  by  about  two  orders  of  magnitude. 
The  technical  details  of  these  rather  cumbersome 
calculations  are  available  in  the  literature  and  will  not 
be  reviewed  here. 

The  results  are  summarized  in  Table  1.  We  show 
estimates  of  all  four  relevant  time  scales  for  the  two 
models  introduced  earlier.  The  “original”  model  [30] 
corresponds  to  nuclear  spins  1  /2  introduced  at  qubits 
in  atoms  without  an  outer  loosely  bound  electron. 
The  “improved”  model  corresponds  to  the  case  when 
the  outer  electron  is  present  and  its  interaction  with 
the  nuclear  spin  and  the  two-dimensional  electron  gas 
dominates  the  dynamics. 

The  data  shown  in  Table  1  were  obtained  assuming 
typical  parameters  for  the  standard  heterojunctions 
utilized  in  quantum-Hall-effect  experiments  today, 
and  qubit  separation  of  65  nm.  Thus,  the  parameter 
values  taken  [28,30]  were  more  appropriate  for  the 
GaAs  system  than  for  Si,  even  though  the  main 
isotopes  of  gallium  and  arsenic  have  nuclear  spin 
3/2  and  cannot  serve  as  spin-zero  hosts.  The  reason 

Table  1 

Time  scales  of  the  qubit  dynamics  for  the  original  [30]  and  improved 
[28]  versions  of  the  nuclear  spin  quantum  computer  with  interac¬ 
tions  mediated  by  the  two-dimensional  electron  gas 


The  original  model 

The  improved  model 

Text 

O(10-5)  sec 

O(10“5)  sec 

Tint 

0(1)  sec 

O(10“2)  sec 

T\ 

O(103)  sec 

0(10)  sec 

t2 

0(10)  sec 

O(l0~')  sec 
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for  using  these  values  has  been  that  experimental 
verification  of  some  of  the  numbers  might  be  possible 
in  the  available  materials  before  cleaner  and  different 
composition  materials  needed  for  quantum  computing 
are  produced. 

Our  estimates,  see  Table  1,  indicate  that  the  quality 
factor  Q  —  10-5  is  not  obtained  for  the  present 
system.  Actually,  no  quantum  computing  proposal 
to  date,  scalable  by  other  criteria,  satisfies  the  10-5 
quality-factor  criterion.  The  values  range  from  1  0~  1  to 
10-2.  The  resolution  could  come  from  development  of 
better  error-correction  algorithms  or  from  improving 
the  physical  system  to  obtain  a  better  quality  factor. 
In  our  estimation  of  the  decoherence  time  scale, 
we  used  parameters  typical  of  a  standard,  “dirty” 
heterostructure  with  large  spatial  fluctuations  of  the 
impurity  potential.  These  heterostructures  have  been 
suitable  for  standard  experiments  because  they  provide 
wider  quantum-Hall  plateaus,  i.e.  ranges  of  magnetic 
field  for  which  all  the  extended  states  of  a  Zeeman 
sublevel  are  filled.  Much  cleaner,  ultra-high  mobility 
structures  can  be  obtained  by  placing  the  ionized 
impurity  layer  at  a  larger  distance  from  the  two- 
dimensional  gas  or  by  injecting  conduction  electrons 
into  the  heterostructure  by  other  means.  Thus,  our 
quantum-computing  proposals  [28,30]  are  unique  not 
only  in  the  large  qubit  separation  allowed  but  also 
in  that  there  is  a  clear  direction  of  exploration  to 
allow  physical,  rather  than  algorithmic,  resolution  of 
the  quality  factor  problem.  This  possibility  should  be 
further  explored  both  experimentally  and  theoretically. 
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Recent  technological  advances  in  electronics  related 
to  spin  polarization  [1,2]  have  boosted  experimental  and 
theoretical  interest  in  quantum  information  science  in 
condensed  matter  systems,  specifically,  in  semiconductor 
heterostructures  at  low  temperatures  and  in  high  magnetic 
fields.  The  solid-state  implementations  of  quantum  in¬ 
formation  devices  seem  to  be  among  the  most  promising 
ones,  due  to  possible  scalability  of  the  elementary  logic 
gates  into  more  complicated  integrated  circuits.  Several 
designs  for  solid  state  and  related  spin-based  quantum 
information  processors  have  been  suggested  [3-8].  Pre¬ 
liminary  experiments,  involving  several  quantum  bits 
(qubits),  have  been  carried  out  or  are  being  contem¬ 
plated  [9,10], 

Our  work  stems  from  the  proposals  that  utilize  nuclear 
or  electronic  spins  as  qubits  for  information  processing 
[3-7].  These  are  natural  choices  for  qubits  because  at  low 
temperatures  spin  states  in  semiconductors  have  relatively 
long  decoherence  times,  sometimes  milliseconds  or  even 
longer  for  electronic  spins,  and  seconds  for  nuclear  spins 
[11-14].  We  propose  a  new  mechanism  for  coupling  be¬ 
tween  two  nuclear-spin  qubits,  combining  aspects  of  two 
models  of  quantum  information  processors,  one  based  on 
nuclear  spins  in  quantum-Hall  effect  systems  [4],  and  an¬ 
other  utilizing  the  nuclear  spins  of  phosphorous  donors  in 
a  silicon  heterostructure  [5]. 

An  appealing  aspect  of  Kane’s  model  [5]  is  a  possibly 
experimentally  feasible  scheme  for  reading  out  the  state  of 
the  quantum  register,  i.e.,  measurement  of  a  nuclear  spin, 
achieved  by  transferring  the  nuclear-spin  polarization  to 
the  electronic  state,  while  the  latter  is  measured  with  the 
use  of  a  single  electron  transistor.  The  model  proposed 
in  [4]  has  a  different  advantage:  Unlike  [5],  the  inter¬ 
action  between  the  nuclear  spins  is  mediated  by  the  two- 
dimensional  (2D)  electron  gas,  and  thus  is  longer  ranged 
due  to  the  highly  correlated  state  of  the  2D  electron  gas 
in  the  quantum-Hall  regime.  This  opens  up  possibilities 
for  experimental  realization  of  such  quantum  information 
processors,  because  large  separation  between  spin  qubits 
means  greater  lithographic  dimensions  in  manufacturing 
the  device.  The  price  paid  is  that  the  coupling  is  weak, 
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and  therefore  the  time  scales  of  the  “gate  function”  can  be 
as  large  as  1  s. 

In  this  work  we  combine  the  two  proposals,  thus  re¬ 
taining  the  measurement  and  control  scheme  proposed  in 
[5,7,9]  and  at  the  same  time  allowing  larger  separations, 
of  the  order  of  100  nm,  between  interacting  qubits.  The 
resulting  system  is  thus  realizable  with  the  present-day 
semiconductor  technologies.  We  propose  a  model  where 
sparsely  positioned  phosphorous  donors  are  imbedded  in 
a  2D  electron  gas  in  the  quantum-Hall  regime.  The  local¬ 
ized  donor  electrons  interact  via  the  delocalized  2D  elec¬ 
trons  and  thus  indirectly  mediate  nuclear-spin  interactions. 
In  3D,  spin  coupling  mechanisms  via  conduction  electrons 
have  been  well  studied  [15].  Here,  we  estimate  the  range 
of  this  induced  nuclear-spin  interaction  for  the  2D  case  and 
find  it  to  be  of  the  order  of  100  nm.  This  is  large  compared 
to  atomic  dimensions,  donor-electron  bound  state  radii, 
and  even  the  electronic  magnetic  length  which  is  typi¬ 
cally  of  the  order  of  10  nm.  We  find  that  this  interaction 
is  also  stronger,  thus  corresponding  to  faster  gate  function 
times,  than  in  [4]. 

We  assume  that  the  coupling  between  the  electronic 
and  the  nuclear  donor  spins  is  given  by  the  Fermi 
contact  interaction,  He-n  =  Acrn  ■  ae.  Here,  A  = 
(87r/3) /Z£,g„/r„|''Lo(0)|2,  where  /i„  and  gn  are  the  nuclear 
magnetron  and  nuclear  g  factor,  respectively,  |’vl/o(0)|2  is 
the  donor-electron  probability  density  at  the  nucleus,  /jlb 
denotes  the  Bohr  magnetron,  and  cr’s  are  Pauli  matrices. 
Coupling  of  the  delocalized  electrons  to  the  nuclear  spin 
is  considerably  weaker  than  that  of  the  localized  donor 
electron.  Therefore,  we  assume  that  the  nuclear  spin 
interacts  with  conduction  electrons  indirectly  via  the 
donor  electron. 

As  a  prototype  system,  we  consider  31P  donors  posi¬ 
tioned  in  Si,  so  all  the  spins  involved  are  The  donor 
electronic  and  nuclear  spins  form  a  four-level  system.  The 
spectrum  of  this  two-spin  system  can  be  obtained  to  O(A) 
with  He~n  treated  as  perturbation.  The  energy  levels  are 
Eo  =  —(yn  +  A)/2  +  A,  E\  =  (y„  -  A)/2  -  A,  E2  = 
(— yn  +  A)/2  —  A,  and  £3  =  ( yn  +  A)/2  +  A,  where 
y„  =  gn/jtnH  is  the  nuclear-spin  splitting.  Here,  H  is  the 
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magnetic  field,  and  the  expression  for  A,  the  electronic 
Zeeman  gap,  will  be  given  shortly.  The  eigenstates  associ¬ 
ated  with  these  energy  levels  are  |0)  =  |  |  e  |  n),  |1)  = 

I  I  e  T  n)  +  (2A/A)  \  \  e  [n),  |2)  =  |  |  e  [  n)  +  (2  A/ 
A)  |  |  e  T  n),  and  |3)  =  |  T  e  T  n),  where  |  |  e  1  n)  rep¬ 
resents  the  electronic  and  nuclear-spin  down  state,  etc. 
Here  we  propose  to  consider  the  states  |0)  and  |1)  as  qubit 
states  of  a  quantum  computer.  By  altering  the  hyperfine 
coupling  constant  A  by  distorting  the  spatial  state  M;o  of 
the  donor  electron  with  an  electrostatic  gate  [5,7],  one 
can  selectively  control  the  state  of  an  individual  qubit  by 
means  of  the  NMR  technique. 

In  order  to  calculate  the  interaction  Hamiltonian  be¬ 
tween  two  qubits,  we  first  consider  the  coupling  between 
the  donor  electron  and  conduction  electrons.  The  ground 
state  of  the  donor  electron  is  bound  (localized)  and  will 
typically  lie  in  the  energy  gap,  several  meV  below  the 
conduction  band  edge.  For  temperatures  of  order  mK, 
electronic  transitions  from  this  localized  state  to  the  con¬ 
duction  band  are  highly  improbable.  The  dominant  in¬ 
teraction  between  the  localized  electron  and  conduction 
electrons  is  their  Coulomb  interaction.  We  are  interested 
only  in  the  exchange  part  of  this  interaction,  i.e.,  the 
spin-dependent  part.  The  spin-independent  part  causes 
screening,  but  it  is  weak  in  2D  [16]  and,  especially  in  the 
presence  of  the  magnetic  field,  cannot  ionize  the  donor. 

In  a  large  magnetic  field,  the  delocalized  2D  electrons 
occupy  highly  degenerate  Landau  energy  levels  [16].  It  is 
convenient  to  introduce  electron  bound  state  creation  and 
annihilation  operators  b}ls  and  bns,  where  n  represents  the 
donor  spatial  state,  and  5  is  the  spin  z  component,  f  or 
].  Let  ahktS,  amkxS  denote  the  creation  and  annihilation 
operators  for  the  delocalized  2D  electrons,  where  m  labels 
the  Landau  level,  while  hkx  is  the  x  momentum  (we  use  the 
asymmetric  gauge).  Then  the  exchange  coupling  between 
the  bound  and  delocalized  electrons  can  be  written  as 


Hex  „  G m.m' .k..k'b\<;amk.s,bn's’Um'k'.s  >  (1) 


where  the  sum  is  over  all  the  indices.  Here,  we  have 
neglected  the  spin-orbit  interaction.  In  what  follows,  we 
will  retain  only  the  lowest  donor-electron  spatial  state,  i.e., 
account  only  for  the  transitions  between  the  two  Zeeman 
levels  of  the  ground  state. 

The  2D  electrons  are  assumed  to  be  in  a  nondissipa- 
tive  quantum-Hall  state  with  filling  factor  v  =  1;  i.e.,  the 
lower  Zeeman  sublevel  of  the  Landau  ground  state  is  com¬ 
pletely  filled  [4].  This  choice  ensures  reduced  decoher¬ 
ence  and  relaxation  effects  [14],  owing  to  the  energy  gap 
in  the  spectrum  of  the  lowest-energy  spin-wave  excitations 
which  are  well  studied  [17,18];  their  spectrum  is  given 


by  £k  =  A  +  Ec[  1  -  /0(€2k2/4)  exp(— €2k2/4)],  where 
Io  is  the  modified  Bessel  function.  Here,  A  =  gpbgH  is 
the  Zeeman  gap,  Ec  =  (tt /7,)]  H(e2  / js  the  character¬ 
istic  Coulomb  energy,  and  g  is  the  effective  g  factor  in  the 
potential  well  that  holds  the  2D  electron  gas,  while  e  is  the 
dielectric  constant  of  the  material,  and  €  =  (fie /e/7) 1/2 
is  the  magnetic  length.  Extension  to  larger  integer  filling 
factors  is  possible  [14,17,18].  One  can  also  introduce  [18] 
normalized  creation  and  annihilation  operators  for  the  spin 
waves,  quadratic  in  electronic  operators, 


ot 

Jk 


+  (kx/2)Aap-(kx/2),h 


(2) 


Here,  Lxy  are  the  transverse  dimensions,  taken  to  infin¬ 
ity  in  the  final  calculation.  The  summation  over  p  is 
taken  in  such  a  way  [18]  that  the  wave  number  subscripts 
are  quantized  in  multiples  of  27 t/Lx.  The  spectrum  of 
these  spin  waves  has  been  experimentally  verified  in  GaAs 
heterostructures  [19]. 

We  will  include  only  these  lowest  excitations  in  the  sum 
(1);  our  goal  is  to  rewrite  (1)  in  terms  of  the  spin-wave 
operators  (2).  The  exchange  coupling  is  thus  truncated  to 

Gm,m',kx,k'x  =  Gkx,k'xS,uo8nifiSmQSmifl,  where 


GkxJi  =  f  ri3Ri</3R2%*(Ri)^o(R2) 

X  U( R,  -  R2)O0\(R2)d>0,*;(Ri),  (3) 


f/(Ri  —  R2)  =  e2/e|Ri  —  R2|  is  the  Coulomb  interac¬ 
tion,  and  'Eo(R)  is  the  donor-electron  ground  state.  The 
states  of  the  conduction  electrons  confined  in  the  2D  well 
are  d>o,jq(R)  =  4>o where  are  the  stan¬ 

dard  2D  Landau  states  [16];  x(z)  describes  the  confine¬ 
ment  of  the  conduction  electron  wave  function  in  the  z 
direction  and  depends  on  the  nature  of  the  confinement 
potential.  Here  and  in  the  following  R  =  (r,  z),  with  R 
and  r  =  (x,y)  being  3D  and  2D  coordinates,  respectively, 
while  z  is  the  direction  perpendicular  to  the  heterostruc¬ 
ture,  in  which  the  applied  magnetic  field  is  pointing. 

With  the  use  of  the  expressions  for  Landau  ground 
state  wave  functions,  </>o,u(r)  =  X 

exp[  — (y  —  f2kx)2 /2(2\.  and  (2),  after  a  lengthy  calcula¬ 
tion,  we  get 

Hex  =  kl  I  e)(l  ,|Sk  +  W*  |  l  e)<T  e\st],  (4) 

L  k 

where  I  I  e)(J  e\  =  b^b\  in  the  appropriate  subspace,  and 

1Vt-«2^U)^-/‘'jR"i3R2^(Rim(R;) 

X  G(Ri  -  R2)Y*(z2)Y(zi)Ck(ri,r2),  (5) 


Ck(ri,r2)  =  expj  - -j-^  [(.ri  -  x2)2  +  (yi  -  y2)2  -  2i(xi  -  x2)(yi  +  y2)]j 


X  exp 


'  A 


(6) 
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Note  that  since  all  the  position  vectors  R,  r  are  measured 
from  the  origin  at  the  donor  atom,  the  quantity  Wk  de¬ 
pends  also  on  the  donor  coordinates.  To  the  leading  order, 
(4)  gives  the  interaction  of  the  donor  electron  spin  with 
excitations  of  the  2D  electron  gas  in  the  v  =  1  integer 
quantum-Hall  state. 

One  can  rewrite  the  interaction  (4) -(6),  with  (4)  multi¬ 
plied  by  the  unit  operator  in  the  nuclear-spin  Hilbert  space, 
in  terms  of  the  eigenstates  of  the  electron-nucleus  system. 
With  the  use  of  the  expressions  derived  earlier  for  these 
eigenstates  in  terms  of  direct  products  of  electronic  and 
nuclear  spin  states,  we  obtain 

Hex  =  }lwk(f|l)<0|  +  |2)(0| 

+  |3)  <1 1  -  y|3)<2|)sk  +  H.c.  (7) 

Now  one  can  calculate  an  effective  Hamiltonian  for  the 
interaction  of  two  qubits.  Since  the  electronic  Zeeman  gap 
is  much  larger  than  the  nuclear  one,  we  can  truncate  the 
Hilbert  space  of  the  combined  electron-nucleus  spins  to 
the  two  lowest  lying  states.  Thus,  we  retain  only  the  |0)  (1 1 
and  conjugate  transitions  in  the  exchange  interaction  (7). 

An  effective  interaction  between  two  qubits  can  be  ob¬ 
tained  within  the  standard  framework  of  second  order  per¬ 
turbation  theory  by  tracing  out  the  states  of  the  spin  waves; 
see  [15,20,21]  for  similar  calculations.  The  result  can  be 
written  as 


This  results  in  control  over  the  effective  coupling  constant 
J  in  (9).  The  precise  effect  of  gates  on  interactions  be¬ 
tween  the  qubits,  as  well  as  on  decoherence  of  their  states, 
should  be  further  studied  in  order  to  establish  the  feasi¬ 
bility  of  the  quantum-computing  approach  proposed  here. 
Most  other  semiconductor  solid-state  quantum-computing 
approaches  [3-7]  utilize  gates. 

Let  us  explicitly  calculate  the  coupling  constant  J  in 

(8)  and  (9).  Because  the  spatial  ground  state  of  the  donor 
is  localized  on  a  scale  smaller  than  the  magnetic  length 
€,  the  overlap  integrand  in  (5)  is  vanishingly  small  for 
lri  —  r2 1  >  i-  At  the  same  time,  for  |k|  >  l/€,  the 
value  of  Ck  decreases  exponentially.  Thus,  Ck  can  be 
simplified  by  neglecting  the  x\  —  x2  and  V|  —  _y2  terms  in 
(6).  Moreover,  for  two  donors  at  separation  larger  than  i, 
we  can  put  (r i  +  r2)/2  —  r;,  with  r7  being  the  location 
of  either  one  of  them.  Then  (5)  can  be  approximated  by 
W]i  j  =  Z(LxLy)  ^H  exp(—  — j-  —  ik  •  r;),  with  Z  = 
(l/27r€2)1/2  /d3R1d3R2^o(Ri)^o(R2)C/(Ri  -R2)  X 
X*(z2)x(zi)- 

Finally,  the  coupling  constant  J  of  the  effective  interac¬ 
tion  (8)  can  be  obtained  by  transforming  the  summation  in 

(9)  to  integration  in  the  limit  Lxy  — > 


J  = 


(M2  |z|2  id 

V  A  /  (2tt)112Ec£2  V  r 


(/■>€), 


(10) 


Hhl  =  7|0, 1 2)  <  1 102 1  +  r|li02)<0il2|.  (8) 

Here,  the  coupling  constant  between  the  two  qubits  is 


/  =  (-)  y 

V  A  /  £k  +  Ex  -  E0 
The  subscripts  1  and  2  in  (8)  and  (9)  label  the  two  donor 
qubits,  while  Wk,i  and  VLki2  are  the  coupling  constants  of 
each  donor  electron  spin  to  spin  waves,  given  by  (5),  and 
£k  is  the  spin-wave  energy. 

The  nuclear-spin  energy  gap  is  much  smaller  than  the 
electronic  spin-wave  excitation  energies.  Therefore,  we 
can  ignore  E\  —  Eq  in  the  denominator  in  (9).  Further¬ 
more,  due  to  the  large  value  of  the  spin-wave  spectral  gap 
at  k  =  0,  £o  =  A,  we  do  not  have  the  “small  denomina¬ 
tor”  problem  encountered  in  other  calculations  of  this  sort, 
e.g.,  [20].  Physically,  this  means  that  the  spin  excitations  in 
the  2D  electron  gas  mediating  the  effective  qubit-qubit  in¬ 
teraction  are  virtual,  and  so  this  interaction  does  not  cause 
appreciable  relaxation  or  decoherence  on  the  gate  function 
time  scale  h/J. 

It  is  important  to  note  that  one  can  construct  a  universal 
CNOT  logic  gate  from  the  controlled  dynamics  governed  by 
Hamiltonians  of  the  form  of  H \_2  and  single  qubit  rotations 
[6].  The  coupling  strength  J  between  the  qubits  can  be 
externally  controlled  by  the  electrostatic  gates  built  above 
the  2D  inversion  layer.  By  applying  gate  voltages,  one  can 
locally  vary  the  density  of  the  2D  electrons,  thus  chang¬ 
ing  coupling  between  the  delocalized  and  donor  electrons. 


wkAwi2 


(9) 


where  d  =  (£'C/2A)1(2£.  A  similar  dependence  of  the 
coupling  on  the  donor  separation  r  was  obtained  in  a 
study  of  nuclear  polarization  diffusion  in  the  quantum-Hall 
regime  [21],  Interaction  (8)  between  the  spins  has  finite 
range  d,  which,  however,  is  very  large  compared  to  the 
effective  Bohr  radius  of  the  donor  ground  state.  Thus, 
the  indirect  exchange  at  large  distances  dominates  the  di¬ 
rect  exchange  interaction  resulting  from  the  overlap  of  the 
two  atomic  wave  functions.  For  magnetic  field  El  =  6  T 
and  e  =  12,  we  get  d  —  65  nm,  which  is  indeed  much 
greater  than  the  characteristic  Bohr  radius  for  a  donor  elec¬ 
tron  in  silicon. 

In  order  to  estimate  J,  we  have  to  evaluate  the  over¬ 
lap  integral  Z.  For  an  order-of-magnitude  estimate,  we 
will  assume  that  x(z)  is  constant  inside  the  well  and 
zero  outside.  Then  Z  —  (2tt)~  [^2{8E)~ 1  /  ri3R]  ri3R2  X 
1L'o(Ri)1Ifo(R2)6?(Ri  “  R2)>  where  8  is  the  width  of  the 
well.  We  put  5—4  nm.  For  M;o(R),  the  donor  ground 
state,  we  choose  a  spherically  symmetric  hydrogenlike 
ground  state  with  the  effective  Bohr  radius  as  —  2  nm. 
This  is,  of  course,  not  the  case  in  a  realistic  situation  [22]. 
The  ground  state  of  the  donor  will  be  influenced  by  the 
band  structure,  by  the  magnetic  field,  and  by  the  confining 
2D  well  potential,  while  the  states  of  the  conducting  elec¬ 
trons  will  be  distorted  by  the  impurity  potentials.  We  are 
not  aware  of  a  thorough  study  of  these  effects  for  our  sys¬ 
tem.  For  the  purposes  of  an  order  of  magnitude  estimate, 
however,  a  spherical  state  should  be  sufficient. 
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Evaluating  the  integral  for  the  Coulomb  potential  U,  we 
obtain  Z  —  (5a2B/\68)Ec.  Assuming  that  the  two  donors 
are  separated  by  the  distance  r  =  100  nm  and  using  the 
value  lA/h  =  58  MHz  from  [4],  we  obtain  the  estimate 
J/h  ~  102  s-1. 

The  clock  speed  of  the  information  processor  just  de¬ 
scribed  appears  to  be  a  fraction  of  kHz  and  should  be 
compared  with  the  time  scales  for  relaxation  and  decoher¬ 
ence.  The  leading  mechanism  for  these  at  low  temperatures 
is  through  interaction  with  impurities.  It  has  been  found 
theoretically  [12,23]  and  confirmed  experimentally  [2]  that 
nuclear-spin  relaxation  in  the  quantum-Hall  regime  is  slow 
and  strongly  dependent  on  the  impurity  potentials;  typi¬ 
cally,  the  relaxation  time  T\  is  of  order  103  s.  In  our  case, 
the  interaction  of  a  qubit  with  the  2D  gas  is  stronger,  and, 
as  a  result,  the  relaxation  is  expected  to  be  faster.  An  es¬ 
timate  from  formulas  in  [12,23]  gives  7j  —  1  s.  There  is, 
however,  another  important  issue — decoherence,  on  time 
scales  To.  Recently,  this  quantity  has  been  calculated  in 
the  same  framework,  that  is,  when  the  interaction  of  the 
conduction  electrons  with  impurities  is  taken  into  account 

[14].  The  results  of  [14]  can  be  adjusted  for  the  present 
case  and  yield  the  estimate  T2  —  10“ 1  s. 

The  existing  quantum  error  correction  protocols  require 
the  quality  factor,  equal  the  ratio  of  the  gate-function  clock 
time  to  decoherence  time,  not  to  exceed  KP5  [24].  Our 
estimates  indicate  that  this  is  not  the  case  for  the  present 
system.  Actually,  no  quantum-computing  proposal  to  date, 
scalable  by  other  criteria,  satisfies  this  10-3  quality-factor 
criterion.  The  values  range  from  10“ 1  to  1(P3.  The  reso¬ 
lution  could  come  from  development  of  better  error- 
correction  algorithms  or  from  improving  the  physical 
system  to  obtain  a  better  quality  factor.  In  our  estimate  of 
the  decoherence  time  scale,  we  used  parameters  typical 
of  a  standard,  “dirty”  heterostructure  with  large  spatial 
fluctuations  of  the  impurity  potential.  These  heterostruc¬ 
tures  have  been  suitable  for  standard  experiments  because 
they  provide  wider  quantum-Hall  plateaus.  Much  cleaner, 
ultrahigh  mobility  structures  can  be  obtained  by  placing 
the  ionized  impurity  layer  at  a  larger  distance  from  the 
2D  gas  or  by  injecting  conduction  electrons  into  the 
heterostructure  by  other  means. 

Thus,  our  present  quantum-computing  proposal  offers 
a  clear  direction  for  exploring  a  physical,  rather  than  al¬ 
gorithmic,  resolution  to  the  quality-factor  problem.  This 
possibility  should  be  further  examined  both  experimentally 
and  theoretically.  Our  new  quantum-computing  paradigm 
suggests  several  interesting  avenues  for  research.  The  ef¬ 
fect  of  gates  on  the  switching  of  qubit  interactions  and  on 
decoherence  requires  further  investigation.  The  first  ex¬ 
perimental  realizations  will  probably  involve  only  a  few 
qubits.  The  interactions  of  these  may  be  significantly 
affected  by  the  geometry,  specifically,  the  edges,  of  the 
heterostructure. 
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Abstract 

We  present  a  new  short-time  approximation  scheme  for  evaluation  of  decoherence. 
At  low  temperatures,  the  approximation  is  argued  to  apply  at  intermediate  times  as 
well.  It  then  provides  a  tractable  approach  complementary  to  Markovian-type  approx¬ 
imations,  and  appropriate  for  evaluation  of  deviations  from  pure  states  in  quantum 
computing  models. 
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1.  Introduction 

Consider  a  microscopic  quantum  system  with  the  Hamiltonian  Hs ■  We  will  refer 
to  the  quantum-computing  single  quantum  bit  (qubit)  or  multi-qubit  paradigm  to  help 
define  the  questions  and  set  up  the  challenges,  in  describing  how  the  system,  S,  interacts 
with  the  surrounding  macroscopic  world.  However,  in  principle  S  can  be  any  quantum 
system. 

Interactions  with  the  surroundings  can  be  quite  different  depending  on  the  setting. 
For  example,  in  quantum  measurement,  which  is  presently  not  fully  understood,  the 
wavefunction  of  the  system  is  probed,  so  part  of  the  process  would  involve  a  strong 
interaction  with  the  measuring  device,  such  that  the  system’s  own  Hamiltonian  plays 
no  role  in  the  process.  However,  in  most  applications,  the  external  interactions  are 
actually  quite  weak.  Furthermore,  the  aim  is  to  minimize  their  effect,  especially  in 
quantum  computing. 

Traditionally,  interactions  with  the  surrounding  world  have  been  modeled  by  the 
modes  of  a  bath,  F>,  with  each  mode  described  by  its  Hamiltonian  Mk,  so  that  the  bath 
of  modes  is  represented  by 


Hb  =  Yjmk-  (1.1) 

K 

The  interaction,  /,  of  the  bath  modes  with  the  system  S,  will  be  modeled  by 

Hi  =  AsPb  =  AsJ2Jk,  (1.2) 

K 

where  A s  is  some  Hermitean  operator  of  S,  coupled  to  the  operator  Pb  of  the  bath. 

The  bath,  or  “heat  bath” ,  can  be  a  collection  of  modes,  such  as  photons,  phonons, 
spins,  excitons,  etc.  For  a  bosonic  bath  of  oscillators,  [1-6],  which  we  use  for  derivation 

-  2  - 
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of  specific  results,  we  take 


Mk  = 


U>Kd'KaK 


(1.3) 


Jk  =  g*KaK  +  gKa]K  .  (1.4) 

Here  we  have  assumed  that  the  energy  of  the  ground  state  is  shifted  to  zero  for  each 
oscillator,  and  we  work  in  units  such  that  h  —  1. 

The  total  Hamiltonian  of  the  system  and  bath  is 


H  =  Hs  +  HB  +  H!  .  (1.5) 

More  generally,  the  interaction,  (1.2),  can  involve  several  system  operators,  each  cou¬ 
pling  differently  to  the  bath  modes,  or  even  to  different  baths.  The  bath  modes,  in  turn, 
can  be  coupled  to  specified  external  objects,  such  as  impurities. 

Let  pit)  represent  the  reduced  density  matrix  of  the  system  at  time  t  >  0,  after 
the  bath  modes  have  been  traced  over.  For  large  times,  the  effect  of  the  environment 
on  a  quantum  system  that  is  not  otherwise  externally  controlled,  is  expected  to  be 
thermalization:  the  density  matrix  should  approach 


p(t  — >  oo) 


exp  j-pHs) 

Tr  s  [exp  (—flHs)]  1 


(1.6) 


where  (3  =  1/kT.  At  all  times,  we  can  consider  the  degree  to  which  the  system  has 
departed  from  coherent  pure-quantum-state  evolution.  This  departure  is  due  to  the 
interactions  and  entanglement  with  the  bath.  We  also  expect  that  the  temperature,  T, 
and  other  external  parameters  that  might  be  needed  to  characterize  the  system’s  density 
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matrix,  are  determined  by  the  properties  of  the  bath,  which  in  turn  might  interact  with 
the  rest  of  the  universe. 

Let  us  introduce  the  eigenstates  of  H$, 

Hs\n)  —  En\n)  ,  (1.7) 

and  have  A E  denote  the  characteristic  energy  gap  values  of  S.  We  also  consider  the 
matrix  elements  of  pit), 


Pmnit )  =  ( m\p(t)\n )  .  (1.8) 

For  large  times,  we  expect  the  diagonal  elements  pnn  to  approach  values  proportional  to 
e~:iE" ,  while  the  off-diagonal  elements,  pm^n ,  to  vanish.  These  properties  are  referred 
to  as  thermalization  and  decoherence  in  the  energy  basis. 

To  establish  these  thermalization  and  decoherence  properties,  several  assumptions 
are  made  regarding  the  system  and  bath  dynamics  [1-11].  At  time  t  —  0,  it  is  usually 
assumed  that  the  bath  modes,  K,  are  thermalized,  i.e. ,  have  density  matrices 


The  density  matrix  R  of  the  system  plus  bath  at  time  t  —  0  is  then  the  direct  product 

*(o)=p(o)n^  (i-io) 

K 

and  the  system  and  bath  modes  are  not  entangled  with  each  other. 

Now,  a  series  of  assumptions  are  made,  e.g.,  the  Markovian  and  secular  approxi¬ 
mations.  The  most  important  is  the  Markovian  approximation,  which,  even  though  it 
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can  be  stated  and  introduced  in  various  ways,  essentially  assumes  that  the  density  ma¬ 
trices  of  the  bath  modes  are  reset  externally  to  the  thermal  ones,  on  time  scales  shorter 
than  any  dynamical  times  of  the  system  interacting  with  the  bath.  This  is  a  natural 
assumption,  because  each  bath  mode  is  coupled  only  weakly  to  the  system,  whereas  it 
is  “monitored”  by  the  rest  of  the  universe  and  kept  at  temperature  T.  In  its  straightfor¬ 
ward  version,  this  amounts  to  using  (1.10)  for  times  t  >  0.  Ultimately,  such  approaches 
aim  at  master  equations  for  the  evolution  of  pmn{t)  at  large  times,  consistent  with  the 
Golden  Rule  and  with  the  expected  thermalization  and  decoherence  properties. 

In  variants  of  these  formalisms,  several  time  scales  are  identified.  One  is  the  inverse 
of  the  upper  cutoff,  Debye  frequency  of  the  bath  modes,  1/ujd-  Another  is  the  thermal 
time  Ti/kT  —  (3  (in  units  of  h  —  1).  The  system  S  has  its  own  characteristic  time, 
1  / A E,  as  well  as  the  system-bath  dynamical  times  of  thermalization  and  decoherence, 
etc.,  Ti,2,...,  corresponding  to  the  “intrinsic”  NMR/ESR  times  Tf,  T2,  etc.  Heuristically, 
bath  modes  of  frequencies  u  comparable  to  A E  are  needed  to  drive  thermalization  and 
decoherence.  Initial  decoherence  can  be  also  mediated  by  the  modes  near  u>  —  0.  At 
low  temperatures,  we  can  assume  that  l/up  <  1/A E  <  (3. 

There  is  evidence  [7,11,12]  that  at  low  temperatures,  the  Markovian-type  and  other 
approximations  used  in  the  derivation  of  equations  for  thermalization  and  decoherence, 
are  only  valid  for  times  larger  than  the  thermal  time  scale  (3.  For  quantum  comput¬ 
ing  applications,  in  solid-state  semiconductor-heterostructure  architectures  [13-19],  we 
expect  temperatures  of  several  tens  of  pK.  The  thermal  time  scale  then  becomes  dan¬ 
gerously  close  to  the  external  single-qubit  control,  Rabi-flip  time  even  for  slower  qubits, 
those  based  on  nuclear  spins.  We  emphasize  that  not  all  the  approximation  schemes 
have  this  problem  [11], 

In  Section  2,  we  offer  additional  comments  on  decoherence  and  quantum  comput¬ 
ing.  Then,  in  Section  3,  we  develop  a  short-time-decoherence  approximation.  In  a 
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discussion  at  the  end  of  Section  3,  we  offer  arguments  that,  at  low  temperatures,  our 
approximation  is  actually  valid  for  intermediate  times,  larger  than  1/u>d,  hopefully  up 
to  times  comparable  or  larger  than  1/A E.  Specific  results  for  the  bosonic  heat  bath  are 
presented  in  Section  4.  Section  5  comments  on  the  case  of  adiabatic  decoherence,  when 
the  short-time  approximation  becomes  exact. 
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2.  Decoherence  and  quantum  computing 

Quantum  computing  architectures  usually  emphasize  systems,  both  the  qubits  and 
the  modes  that  couple  them  (and  at  the  same  time  act  as  a  bath  mediating  unwanted 
coupling  to  the  rest  of  the  universe),  that  have  large  spectral  gaps.  It  is  believed  that, 
especially  at  low  temperatures,  spectral  gaps  slow  down  relaxation  processes.  Therefore, 
quantum  computing  architectures  usually  assume  [13-19]  qubits  in  quantum  dots,  or  in 
atoms,  or  subject  to  large  magnetic  fields,  and  coupled  by  highly  nondissipative  quantum 
media  [14,19]. 

The  spectral  gaps  are  expected  to  slow  down  exponentially,  by  the  Boltzmann 
factor,  the  processes  of  thermalization,  involving  energy  exchange.  Off-shell  virtual 
exchanges,  will  be  also  slowed  down,  but  less  profoundly.  The  latter  processes  contribute 
to  decoherence.  Therefore,  at  low  temperatures,  we  might  expect  separation  of  time 
scales  of  the  initial  decoherence  vs.  later-stage  thermalization  and  further  decoherence. 
The  latter  two  processes  are  described  by  the  traditional  NMR/ESR  intrinsic  Tj  and 
X2,  respectively. 

Since  only  thermalization  is  clearly  associated  with  the  energy  eigenbasis,  one  can 
also  ask  whether  the  energy  basis  is  the  appropriate  one  to  describe  decoherence  for 
short  and  intermediate  times,  before  the  thermalizing  processes,  that  also  further  drive 
decoherence,  take  over.  The  issue  of  the  appropriate  basis  for  studying  decoherence,  has 
also  come  up  in  models  of  quantum  measurement.  It  has  been  argued  [20-24]  that  the 
eigenbasis  of  the  interaction  operator,  Ag,  may  be  more  appropriate  for  intermediate 
times  than  the  energy  eigenbasis. 

Yet  another  aspect  of  decoherence  in  quantum  computing,  involves  the  observation 
that  we  really  want  to  retain  a  pure  state  in  the  quantum  computation  process  [25- 
30].  Decay  of  off-diagonal  matrix  elements,  in  whatever  basis,  might  not  be  the  best 
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measure  of  deviations  from  the  pure-state  density  matrix.  For  instance,  the  deviation 
of  Tr  s  [p2  (£)]  from  1,  may  be  more  appropriate.  Therefore,  it  is  desirable  to  have 
basis-independent  expressions  for  the  reduced  density  operator  p(t). 

Recently,  several  groups  have  reported  [12,19,24,31-41]  results  for  spin  decoherence 
in  solid  state  systems  appropriate  for  quantum  computing  architectures.  Some  of  these 
works  have  not  invoked  the  full  battery  of  the  traditional  approximations,  Markovian 
and  secular,  etc.,  or  have  utilized  the  spectral  gap  of  the  bath  modes,  to  achieve  better 
reliability  of  the  short-time  results.  In  [41],  interaction  of  the  spin-exciton  bath  modes 
with  impurities  was  accounted  for,  as  the  main  mechanism  of  decoherence.  In  the 
present  work,  we  limit  ourselves  to  the  bath  modes  only  interacting  with  the  system. 
Experimental  efforts  are  picking  up  momentum,  with  the  first  limited  results  available 
[42,43]  by  traditional  NMR/ESR  techniques,  with  the  quantum-computing  emphasis. 

An  approach,  termed  adiabatic  decoherence,  have  been  developed  by  us  [24],  ex¬ 
panding  the  earlier  works  [12,31-33],  with  the  goal  of  avoiding  the  ambiguity  of  the  basis 
selection  and  achieving  exact  solvability.  The  price  paid  was  the  assumption  that  H s 
is  conserved  (a  particular  version  of  the  quantum  nondemolition  processes),  which  is 
equivalent  to  requiring  that 

[Hs,H]  =  [Hs,  As]  =  0  (adiabatic  case)  .  (2-1) 

This  makes  the  eigenbasis  of  Hs  and  A g  the  same,  but  precludes  energy  relaxation, 
thus  artificially  leaving  only  energy-conserving  relaxation  pathways  that  contribute  to 
decoherence.  We  will  comment  on  the  results  of  this  approach  in  Section  5. 

Most  of  the  results  referred  to  earlier,  have  involved  approximations  of  one  sort  or 
another.  The  most  popular  and  widely  used  approximation  has  been  the  second-order 
perturbative  expansion  in  the  interaction  strength,  Hi,  though  some  nonperturbative 
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results  have  also  been  reported.  In  Section  3,  we  describe  a  novel  approximation  scheme 
[44]  that  is  valid  for  short  times.  It  has  several  advantages,  such  as  becoming  exact  in 
the  adiabatic  case,  allowing  derivation  of  several  explicit  results,  and,  at  least  in  princi¬ 
ple,  permitting  derivation  of  higher-order  approximations.  Certain  models  of  quantum 
measurement  evaluate  decoherence  by  effectively  setting  II  s  =  0.  Our  approximation 
then  becomes  exact,  and  our  results  are  consistent  with  these  studies  [45,46]. 

Our  formulation  in  Section  3,  will  be  quite  general,  and  we  will  not  use  the  specific 
bath  or  thermalization  assumptions.  However,  we  do  utilize  the  factorization  property 
(1.10)  at  time  t  —  0.  Thus,  we  do  have  to  assume  that,  at  least  initially,  the  system  and 
the  bath  modes  are  not  entangled.  In  fact,  the  present  formulation  also  relies  on  that 
the  Hamiltonians  at  hand  are  all  time-independent.  Therefore,  we  have  excluded  the 
possibility  of  controlled  dynamics,  in  the  quantum  computing  sense,  when  gate  functions 
are  accomplished  by  external  couplings  to  individual  qubits  and  by  external  control  of 
their  pairwise  interactions.  Our  formulation,  therefore,  applies  to  “idling”  qubits  or 
systems  of  (possibly  interacting)  qubits.  It  is  reasonable  to  assume  that  a  lower  limit 
on  decoherence  rate  can  be  evaluated  in  such  an  idling  state,  even  though  for  quantum 
error  correction,  qubits  otherwise  idling,  might  be  frequently  probed  (measured)  and 
entangled  with  ancillary  qubits  [25-30]. 

The  t  =  0  factorization  assumption  (1.10),  shared  by  all  the  recent  spin-decoherence 
studies,  then  represents  the  expectation  that  external  control  by  short-duration  but  large 
externally  applied  potentials,  measurement,  etc.,  will  “reset”  the  qubits,  disentangling 
them  from  the  environment  modes  to  which  the  affected  qubits  are  only  weakly  coupled. 
Thus,  we  assert  that  it  is  the  qubit  system  that  gets  approximately  reset  and  disentan¬ 
gled  from  the  bath  towards  time  t  —  0,  rather  than  the  bath  is  thermalized  by  the  rest 
of  the  universe,  as  assumed  in  Markovian  approximation  schemes. 
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3.  Short-time  decoherence 

In  addition  to  the  energy  basis,  (1.7),  we  also  define  the  eigenstates  of  the  interaction 
operator  As,  by 


AsIt)  =  A7|y)  ,  (3.1) 

where  the  Greek  index  labels  the  eigenstates  of  As,  with  eigenvalues  A7,  while  the 
Roman  indices  will  be  used  for  the  energy  basis,  and,  when  capitalized,  for  the  bath 
modes,  (1.2)- (1.4). 

The  time  dependence  of  the  density  matrix  R{t )  of  the  system  and  bath,  is  formally 
given  by 


R(t)  =  eri{Ha+HB+H,)tR( 0)  e‘(Hs+HB+H,)t 


(3.2) 


We  will  utilize  the  following  approximate  relation  for  the  exponential  factors,  as  our 
short-time  approximation, 


^(Hs  +  Hs+H^t+Oit3)  _  eiHst/2  ei{HB+H!)t  eiHst/2 


(3.3) 


This  relation  has  the  following  appealing  properties.  It  becomes  exact  for  the  adiabatic 
case,  (2.1).  Furthermore,  if  we  use  the  right-hand  side  and  its  inverse  to  replace  e±lHt , 
then  we  are  imposing  three  time-evolution-type  transformations  on  i?(0).  Therefore, 
the  approximate  expression  for  R(t)  will  have  all  the  desired  properties  of  a  density 
operator.  Finally,  extensions  to  higher-order  approximations  in  powers  of  t  are  possible, 
by  using  relations  derived  in  [47],  where  various  expressions  valid  to  0(t 4)  and  0(t 5) 
were  considered. 
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Our  goal  is  to  evaluate  the  resulting  approximation  to  the  matrix  element, 


pmn(t)  =  Tr  B(m\e-‘Hst,‘‘  e-‘(HB+H,>t  e~iH3t/2R{0)  e‘Hst/2  eMHB+Hl)t  e‘Hst/2\n}  . 

(3.4) 

First,  we  apply  the  operators  Hs  in  the  outer  exponentials,  acting  to  the  left  on  (m\  , 
and  to  the  right  on  | n),  replacing  II s  by,  respectively,  Ern  and  En.  We  then  note  that 
the  second  exponential  operator  in  (3.4)  contains  Ag,  see  (1.2).  Therefore,  we  insert  the 
decomposition  of  the  unit  operator  in  the  system  space,  in  terms  of  the  eigenbasis  of  Ag, 
before  the  second  exponential,  and  one  in  terms  of  the  eigenbasis  of  iJg  after  it.  This 
allows  us  to  apply  A g  in  the  second  exponential  and  also  II  s  in  the  third  exponential. 
The  same  substitution  is  carried  out  on  the  other  side  of  i?(0),  with  the  result 


Pmn(t )  —  ^  '  Tt  b 

7  p  q  S 


'™\l)(l\p)e 


o) 


x  e^HB+xsPB^>t(ci\^)(^\n)eiE'nt^2 


(3.5) 


K 


The  next  step  is  to  collect  all  the  terms,  and  also  identify  that  the  trace  over  the 
bath  can  be  now  carried  out  for  each  mode  separately.  We  use  (1.1)-(1.2)  to  write 


Pmn(t)  =  {el{Eq+En  Ep  Em)t/2(m\-i)(~i\p)  ppq(tt)  <#)(%) 

7  p  q  S 


X 


n 


Tr 


K 


(3.6) 


K 

While  this  expression  looks  formidable,  it  actually  allows  rather  straightforward  calcula¬ 
tions  in  some  cases.  Specifically,  the  simplest  quantum-computing  applications  involve 
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two-state  systems.  Therefore,  the  sums  in  (3.6)  are  over  two  terms  each.  The  calcula¬ 
tions  involving  the  overlap  Dirac  brackets  between  the  eigenstates  of  H$  (labeled  by  m, 
n,  p  and  q)  and  of  A s  (labeled  by  7  and  5),  as  well  as  the  energy-basis  matrix  elements 
of  p(0),  cf.  (1.8),  involve  at  most  diagonalization  of  two-by-two  Hermitean  matrices. 
Of  course,  the  approximation  (3.6)  can  be  used  for  evaluation  of  short-time  density 
matrices  for  systems  more  general  than  two-state. 

The  challenging  part  of  the  calculation  involves  the  trace  over  each  mode  of  the 
bath.  Since  these  modes  have  identical  structure,  e.g.,  (1.3)-(1.4)  for  the  bosonic  bath 
case,  but  with  AT-dependent  coupling  constants,  the  calculation  needs  only  be  done 
once,  in  the  space  of  one  mode.  Furthermore,  results  for  the  bath  models  ordinarily 
used,  such  as  the  bosonic  and  spin  baths,  are  either  already  available  in  the  literature 
or  can  be  calculated  without  much  difficulty.  For  the  thermalized  initial  bath-mode 
density  matrix  9k,  we  give  the  exact  bosonic-model  expression  in  the  next  section. 

In  the  remainder  of  this  section,  we  first  further  analyze  the  trace  over  one  bath 
mode  entering  (3.6).  We  then  comment  on  the  limits  of  validity  of  the  present  approx¬ 
imation. 

In  an  obvious  shorthand  notation,  we  write  the  single-mode  trace  in  (3.6)  as 


Tr 


e 


Tr 


0ei(M+5J)te 


(3.7) 


Now,  to  the  same  order  of  approximation  as  used  in  (3.3),  we  can  write 


gi(M+8J)t+0(t3)  _  ^ iMt/2  giSJt  iMt/2 


The  resulting  approximation  for  the  trace  (3.7)  reads 


(3.8) 
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Tr 


iMt/2  q 


ei((5-7)Jt 


(3.9) 


which  illustrates  that,  within  this  approximation,  the  product  of  traces  in  (3.6)  is  a 
function  of  the  difference  A7  —  As-  In  fact,  this  product  is  exactly  1  for  A7  =  A s  and,  in 
most  applications,  the  following  form  is  likely  to  emerge, 


Yl  Tr K [...]  =  e-const(A7-A.)2t2+o(i3)  ^  (3.10) 

K 

though  we  caution  the  reader  that  (3.10)  is  somewhat  speculative  and  suggested  by  the 
exact  result  for  the  bosonic  heat  bath,  reported  in  the  next  section. 

Finally,  we  point  out  that  in  most  cases  of  interest,  the  initial  single-mode  density 
matrix  6  will  commute  with  the  bath-mode  energy  operator  M.  In  fact,  the  thermalized 
6  is  a  function  of  M.  Therefore,  (3.9)  can  be  further  simplified  to 


Tr 


Q  e i{S-j)Jt 


(3.11) 


However,  let  us  emphasize  that  the  approximate  relations  (3.9)-(3.11)  are  likely  of 
value  only  as  far  as  they  help  to  derive  basis-independent  (operator)  approximations  to 
pit),  by  a  technique  illustrated  in  the  next  section.  Indeed,  for  most  bath  models  it  is 
advisable  to  calculate  the  single- mode  trace  exactly  first,  according  to  (3.6),  and  then 
attempt  various  approximations. 

The  latter  statement  reflects  our  expectation  that  the  approximation  developed 
here  is  valid,  for  low  temperatures,  not  only  for  short  times,  defined  by  t  <  1/u>d, 
but  also  for  intermediate  times,  exceeding  1/u>d-  This  is  suggested  by  the  result  of 
an  illustrative  calculation  in  the  next  section,  but  mainly  by  the  fact  that  (3.11)  only 
includes  the  bath-mode  energy  scales  via  6 ,  and,  therefore,  at  low  temperatures,  is 
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dominated  by  the  lowest  bath-mode  excitations,  and  is  not  sensitive  to  frequencies 
of  order  lod-  Thus,  we  expect  our  approximation  to  be  applicable  complementary  to 
the  Markovian-type  approximations  and  definitely  break  down  in  the  regime  of  fully 
developed  thermalization,  for  t  >  ()(j3).  Additional  supporting  observations  are  offered 
in  Section  5,  when  we  consider  the  adiabatic  case  (2.1). 
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4.  The  bosonic  heat  bath 

In  this  section,  we  consider  the  bosonic  heat  bath  [6],  see  (1.3)- (1.4),  in  the  initially 
thermalized  state, 


9k  =  g~^Mk I  Tr  K  (e_/3Mif)  =  (l  -  e~l3u}K)  .  (4.1) 

The  product  of  the  single- mode  traces  in  (3.6),  is  then  available  in  the  literature 
[12,24,31], 


pmn(t)  =  Y  |el(^+s"  E p  Em)t/2(m\^)(-f\p)(q\S)(6\n)ppq{0) 
7  p  q  5 


x  exp 


-E 


\9k\ 


2  r 


K 


id 


K 


o  ,\  n2  .  2  UKt  @U)K 

2  —  As )  sm  — —  cotli  — - — 


+  i  (A^  -  X2S)  (smuKt  -  uKt)  ) 


(4.2) 


The  last  term  in  the  exponent,  linear  in  t,  is  usually  viewed  as  “renormalization”  of  the 
system  energy  levels  due  to  its  interaction  with  the  bath  modes.  It  can  be  removed  by 
adding  the  term, 


Hr  =  A|  Y  \9k\2/uk  ,  (4.3) 

K 

to  the  total  Hamiltonian.  However,  the  usefulness  of  this  identification  for  short  times 
is  not  clear,  and  we  will  not  use  it.  One  can  check  that,  unmodified ,  (4.2)  is  consistent 
with  the  expectation  (3.10). 

Let  us  now  define  two  non- negative  real  spectral  sums,  B(t)  and  C(t).  over  the 
bath  modes, 
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(4.4) 


(4.5) 


When  converted  to  integrals  over  the  bath  mode  frequencies,  with  the  cutoff  at  up, 
these  sums  have  been  discussed  extensively  in  the  literature  [6,12,31],  for  several  choices 
of  the  bath  mode  density  of  states  and  coupling  strength  g  as  functions  of  the  mode 
frequency. 

The  final  expression  is, 


x  exp  -jz?2(l)  (A7  -  A5)2  -  iC(t)  (A2  -  A^) 


(4.6) 


When  the  spectral  functions  are  expanded  in  powers  of  t,  this  result  confirms  all  the 
conclusions  and  conjectures  discussed  in  Section  3,  in  connection  with  relations  (3.9)- 


(3.11). 


Let  us  now  turn  to  the  derivation  of  the  basis-independent  representation  for  pit), 
by  utilizing  the  integral  identity 


(4.7) 


Exponential  factors  in  (4.6)  can  then  be  reproduced  by  applying  operators  on  the  wave- 
functions  entering  the  overlap  Dirac  brackets,  with  the  result 
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iHst/2  ei[yB(t)As  —  C(t)Ag]  g  —  iHst/2  p(Q\  g* Hst/2  e~i[yB(t)As  —  C(t)A2]  g iHst/2 

(4.8) 


Within  the  0(t 2)  approximation  (3.3),  given  that  B  and  C  are  of  order  linear  or 
higher  in  t,  we  can  combine  the  exponential  operators  to  get  an  alternative  approxima¬ 
tion, 


V*  pit)  = 

though  (4.6)  and  (4.8)  are  in  fact  easier  to  handle  in  actual  calculations. 

As  an  application,  let  us  consider  the  case  of  Hs  proportional  to  the  Pauli  matrix 
a2,  e.g.,  a  spin-1/2  particle  in  magnetic  field,  and  A s  —  crx,  with  the  proportionality 
constant  in  the  latter  relation  absorbed  in  the  definition  of  the  coupling  constants  g k 
in  (1.4).  Let  us  study  the  deviation  of  the  state  of  a  spin- 1/2  qubit,  initially  in  the 
energy  eigenstate  |  j )  or  |  j. ),  from  pure  state,  by  calculating  Tr  $  [p2{t)\  according  to 
(4.8).  We  note  that  for  a  two-by-two  density  matrix,  this  trace  can  vary  from  1  for  pure 
quantum  states  to  the  lowest  value  of  1/2  for  maximally  mixed  states. 

A  straightforward  calculation  with  p(0)  =  |T)(t|or|j)(j|,  yields 


dye 


-y2  - i[tHs-yB(t)As+C(t)A |] 


p{  0)e 


i[tHs-yB(t)As+C(t)A2\ 


(4.9) 


Tr  s[p2(t)] 


1 

2 


1  +  e 


2  B2(t) 


(4.10) 


As  the  time  increases,  the  function  B2{t )  grows  monotonically  from  zero  [6,12,24,31]. 
Specifically,  for  Ohmic  dissipation,  B2(t )  increases  quadratically  for  short  times  t  < 
0(1/ud),  then  logarithmically  for  0(1/ud)  <  t  <  0(K/kT),  and  linearly  for  t  > 
0{J%/kT).  (For  other  bath  models,  it  need  not  diverge  to  infinity  at  large  times.)  This 
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calculation  thus  illustrates  the  fact  that  the  present  approximation  can  yield  reasonable 
results  for  short  and  even  intermediate  times. 

Both  approximations,  (4.8)-(4.9),  make  the  deviation  from  a  pure  state  p(0)  = 
IV’oXV’ol  apparent:  p(t  >  0)  is  obviously  a  mixture  (integral  over  y )  of  pure-state  pro¬ 
jectors  \ijj(y,t)){^(y,t) |,  where,  for  instance  for  (4.9), 

tp(y,  t )  =  e-*Hs-yB(t)As+C(t)A%]  ^  ?  (4.H) 

with  a  somewhat  different  expression  for  (4.8). 
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5.  The  adiabatic  case 

Relation  (2.1)  corresponds  to  the  system’s  energy  conservation.  Therefore,  energy 
flow  in  and  out  of  the  system  is  not  possible,  and  normal  thermalization  mechanisms  are 
blocked.  The  fact  that  our  approximation  becomes  exact  in  this  case,  provides  support 
to  the  expectation  that,  at  low  temperatures,  it  is  generally  valid  beyond  the  cutoff  time 
scale  1/Vd,  providing  a  reasonable  evaluation  of  decoherence  and  deviation  from  a  pure 
state,  as  exemplified  by  the  calculation  yielding  (4.10),  in  Section  4. 

With  (2.1),  we  can  select  a  common  eigenbasis  for  Hs  and  As-  Then  the  distinction 
between  the  lower-case  Roman  and  Greek  indices  in  (3.6)  becomes  irrelevant,  and  the 
sums  can  all  be  evaluated  to  yield 


pmn  (t) 


—  g  *(-®n 


Em)tprnn{ 0)  Yl  Tl'K 
K 


e  —  g*(Afjf+X„  Jif)t 


(5.1) 


This  expression  was  discussed  in  detail  in  our  work  on  adiabatic  decoherence  [24],  Specif¬ 
ically,  for  the  initially  thermalized  bosonic  heat  bath  case,  we  have,  for  the  absolute 
values  of  the  density  matrix  elements, 


\pmn(t)\  =  |pmn(0)|  e-B2(mm-Xn)2/4  (5.2) 

The  decay  of  the  off-diagonal  matrix  elements  thus  depends  of  the  properties  of  the 
spectral  function  B2(t)  as  the  time  increases.  Such  explicit  results  [12,24,31-33]  illustrate 
that  for  true  irreversibility,  the  number  of  bath  modes  must  be  infinite,  with  the  spectral 
function  evaluated  in  the  continuum  limit. 

In  summary,  we  have  derived  short-time  approximations  for  the  density  matrix  and 
its  energy-basis  matrix  elements.  Our  expressions  are  quite  easy  to  work  with,  because 
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for  few-qubit  systems  they  only  involve  manipulation  of  finite-dimensional  matrices,  and 
they  will  be  useful  in  estimating  decoherence  and  deviation  from  pure  states  in  quantum 
computing  models,  including  results  for  low  temperatures. 
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